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Abstract 

We reexamine the theoretical limit on the axion coupling to two photons, Qa-f-f, in 
the light of a new solar model (seismic solar model) and the latest solar neutrino 
observations, the Super-Kamiokande and the Sudbury Neutrino Observatory (SNO). 
From the comparison of the theoretically expected and the measured neutrino fluxes, 
we set a limit ga-yy < 4.0 x 10“^^ GeV~^. This limit based on a new procedure is 
about a factor of 3 improvement over the previous theoretical limit and a more severe 
limit than the solar axion experiments. Therefore this limit is the most stringent 
limit on solar axions. 
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1 Introduction 


The axion is a light pseudoscalar particle introduced to solve the strong CP 
problem [1,2] and one of the most likely candidates for the dark matter. There¬ 
fore the axion is an extremely important particle for particle physics and as¬ 
tronomy. 

In the interior of the Sun, blackbody photons can convert into axions in 
the fluctuating Coulomb helds of the charged particles in the plasma, 7 + 
{e~,Ze) {e~,Ze) + a, and this reaction is called as the Primakoff pro¬ 

cess. Large number of current and proposed experiments are trying to detect 
thus-created solar axions [3-6]. 

The axion interacts so weakly with other particles that it escapes freely from 
the Sun once it is produced. Therefore it functions as an energy-loss mecha¬ 
nism. Schlattl et ah [7] calculated evolutionary solar models taking account of 
this axionic energy-loss. Comparing those models with constraints from helio- 
seismology and neutrino oscillations, they derived a solar limit on the axion 
coupling to two photons < 10.0 x 10“^° GeV“^ 

There are, however, some complaints about their approach as follows: 

• Evolutionary solar models are constructed in the framework of stellar evo¬ 
lution, in which various assumptions about the past history of the Sun have 
to be adopted. 

• Constraints from helioseismology are not effectively utilized. 

Therefore we reexamine the calculation by using seismic solar models, which 
are free from evolutionary assumptions and utilize helioseismic constraints to 
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the full. 


2 Seismic solar model 

Helioseismology has been successful in determining precisely the sound-speed 
prohle, c(r), and the density prohle, p(r), in the Sun and the depth of the con¬ 
vection zone, Tconv, based on precise observations. By imposing the constraints 
of Cobs(^"), Pobs(^"), and Tcom, wc can solve the basic equations governing the ra¬ 
diative core of the Sun (the continuity equation, the hydrostatic equation, the 
energy equation, and the energy transfer equation) directly without following 
the evolutionary history of the Sun. We call this model as the seismic solar 
model (hereafter SeiSM) and SeiSM is described in detail in Ref. [8]. SeiSM 
has advantages over the conventional solar models, evolutionary models, as 
follows: 

• We can construct a model of the present-day Sun and evaluate the theoreti¬ 
cally expected neutrino fluxes without assuming and following the evolution¬ 
ary history of the Sun, which cannot be justified directly by observations. 

• We can construct a model without worrying about the treatment of convec¬ 
tion [9], which are not well-described theoretically. 

• SeiSM is faithfully consistent with almost all observations (cobs(’"), Pobs(^), 
Tconv, luminosity, and the mass ratio of heavy elements to hydrogen at 
the surface, which is determined spectroscopically) except for the neutrino 
fluxes, while evolutionary models are not necessarily so, as shown later in 
Fig. 1 and 2. 
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3 Seismic solar models with axionic energy-loss 


The axionic energy-loss rate by the Primakoff effect can be written in the form 
^axion = 0.892 X 10"^ qIq T/ P2 ^ ^rg g“^ s"\ (1) 


F{k^) 








( 2 ) 


~ 8.28p2Tf3 (3 + X), 


(3) 


where pio = Pa.y.y/10“^° GeV“^, Tj = T/l(f K, p2 = p/10^ g cm“^, and X 
is the mass fraction of hydrogen [7,10,11]- We calculate a series of SeiSMs 
with varying Pa 775 which determines the amount of axionic energy-loss. The 
model’s properties are summarized in Table 1. We show the axion luminosity. 
La, the temperature, T, and the mass fraction of heavy elements, Z, at the 
core, the mass fraction of helium, Y, at the surface, and the theoretically 
expected ®B-neutrino flux and neutrino capture rates for the chlorine and the 
gallium experiments. In this procedure, SeiSMs are constructed so as to be 
always consistent with Cobs(^") and Pobs(^") (Fig. 1 and 2). In addition, Wurf of 
thus-constructed SeiSMs are always in satisfactory agreement with the values 
determined directly from helioseismic inversions as demonstrated in Ref. [8]. 

We do not try to explain how SeiSM with axionic energy-loss can be realized in 
the evolutionary process, and such an investigation is beyond our scope. Our 
purpose is to clarify what can be derived if we construct self-consistent solar 
models with axionic energy-loss, which are faithfully consistent with almost 
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all observations except for the neutrino fluxes and free of severe subjective 
restrictions on the evolutionary history. In our opinion, assumptions in the 
standard evolutionary solar models (e.g. initial uniform chemical composition, 
no mass loss, and no accretion) are subjective restriction, which cannot be 
justihed directly by observations. 

An increase in the axionic energy-loss must be compensated by an increase in 
the nuclear energy generation, which leads to an increase in the theoretically 
expected neutrino fluxes. An increase in nuclear energy generation is realized 
with increases in T and p near the center of the Sun. However, p of SeiSM is 
constrained by observations, so T should increase particularly. A higher Tcore 
means a steeper temperature gradient, which is proportional to the opacity. 
A higher opacity is realized with a higher Z. All these trends can be seen in 
Table 1, Fig. 3, and later in Fig. 4. Although Z-profiles shown in Fig. 3 may 
appear rough, it is sufficient to calculate the Z-pro£le of SeiSM to the third 
decimal place for discussing most of the properties of SeiSM other than the 
Z-prohle, itself [8]. 

The density constraint and a non-uniform Z-prohle make Tcore of SeiSM very 
sensitive to Qa'y'y- Because the nuclear reaction rates are mainly controlled by 
Tcore, the theoretically expected neutrino fluxes of SeiSM are more sensitive to 
than those of evolutionary solar models with a uniform Z-pro£le [7]. We 
found that if we construct seismic solar models with a uniform Z-profile [14], 
which are not necessarily consistent with Pobs(’"), by imposing the constraints 
of Cobs(^") and rconv, responses of the theoretically expected neutrino fluxes of 
such models to pa'j'y are similar to those of evolutionary solar models with a 
uniform Z-pro£le [7]. 
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4 Solar neutrinos 


Recently the Sudbury Neutrino Observatory (SNO) has announced the hrst 
results on solar neutrinos [15]. SNO detected the ®B-neutrino via the charged 
current (CC) reaction on deuterium and by the elastic scattering (ES) of elec¬ 
trons. The CC reaction is sensitive exclusively to z/e, while the ES reaction, 
the same reaction as the Super-Kamiokande (SK), is sensitive to all active 
neutrino flavors {ue, i^t), but with weak sensitivity to and z/,-. If there 
are flavor transformations between active neutrino flavors (z/e, z/^, z/,-), compar¬ 
ison of the neutrino flux deduced from the ES reaction assuming no neutrino 
oscillations, ^^^(z/j.), to that measured by the CC reaction, ^'"‘"(z/e), can pro¬ 
vide the flux of non-electron flavor active neutrinos, t)- By comparing 
to the Super-Kamiokande’s precise value of 0^®(z/j,) [16], the total flux 
of the active ®B-neutrino, 0(z/e, r), is determined to be 5.44±0.99xl0®cm“^ 

s“^ [15]. If there are also oscillations to sterile neutrinos, the total flux may 
be slightly larger. 

The neutrino flux derived in this way is consistent with the theoretically ex¬ 
pected neutrino flux of SeiSM with ga-y-y = 0. The relation of the observa- 
tionally determined and the theoretically expected neutrino fluxes are shown 
in Fig. 4. We estimate the uncertainty of the theoretically expected neutrino 
flux by the method used in Ref. [8]. An increase in Qa-y-y leads to an increase in 
the theoretically expected neutrino flux of SeiSM as explained in the previous 
section. Because a higher Qayy than a critical value makes the model’s neutrino 
flux too high to be consistent with that determined from SNO and SK, we 
can limit ga^yy- From the comparison of these neutrino fluxes, we set a limit 
ga'y'y < 4.0 X 10“^° GeV“^. 
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5 Conclusion 


We calculate a series of seismic solar models (SeiSM) with varying the ax- 
ion couplings to two photons, Qa'y'y, which determines the amount of axionic 
energy-loss. An increase in axionic energy-loss leads to increases in the nu¬ 
clear reaction rates and the expected neutrino fluxes. The theoretically ex¬ 
pected ®B-neutrino flux of SeiSM should be identical with the total flux of 
the active neutrino flux, 0(z/e, r), which is determined by comparing the 

measured neutrino flux of the Sudbury Neutrino Observatory and that of the 
Super-Kamiokande. From the comparison of these neutrino fluxes, we set a 
limit Qa'y'y < 4.0 X GeV“h This limit based on a new procedure is about 
a factor of 3 improvement over the previous theoretical limit [7] and a more 
severe limit than the solar axion experiments [3,4]. Therefore this limit is the 
most stringent limit on solar axions. 
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Fig. 1. Relative differences in the square of the sound-speed between seismic solar 
models with axionic energy-loss and the helioseismically determined profile [12], 
c^(obs). For a comparison, the latest evolutionary model, BPOO [13], is also shown. 



Fig. 2. Relative differences in the density. The line styles and references are the 
same as given in Fig. 1. 
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Fig. 3. Z-profiles of seismic solar models with different qiq = 5(a77/10 GeV ^ 



Fig. 4. The theoretically expected ®B-neutrino fluxes of seismic solar models with 
different 5 a 77 - For a comparison, the observationally determined active neutrino 
flux is also shown by the horizontal solid line with the ±1 a error band (the dotted 
lines). Vertical lines are upper limits for this work, the previous theoretical 

limit [7], and the most stringent limit from the solar axion experiments [4]. 
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Table 1 

Seismic solar models with axionic energy-loss 


5a77 (10-^° GeV-1 

) 0.0 

3.0 

4.0 

4.5 

Aa/ Lq 

0.000 

0.019 

0.037 

0.050 

Tcore (10^ K) 

1.58 

1.61 

1.63 

1.65 

^core 

0.022 

0.027 

0.031 

0.035 

Isurf 

0.240 

0.240 

0.240 

0.240 

Neutrino detection 

rates 




(10® cm“^ s“^ ) 

5.9 

8.3 

10.8 

13.4 

Cl (SNU*) 

8.7 

11.8 

14.9 

18.2 

Ga (SNU) 

132 

147 

161 

175 


* A SNU is defined to be 10 interactions s ^ per target atom. 
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